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We analyze the quantum entanglement at the equilibrium in a class of exactly solvable one- 
dimensional spin models at finite temperatures and identify a region where the quantum fluctuations 
determine the behavior of the system. We probe the response of the system in this region by studying 
the spin dynamics after projective measurement of one local spin which leads to the appearance of 
the "decoherence wave" . We investigate time-dependent spin correlation functions, the entanglement 
dynamics, and the fidelity of the quantum information transfer after the measurement. 
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I. INTRODUCTION 

Collective behavior in many-body quantum systems is 
associated with the development of classical correlations, 
as well as of the correlations which cannot be accounted 
for in terms of classical physics, namely, entanglement. 
The entanglement represents in essence the impossibility 
of giving a local description of a many-body quantum 
state. Experimental tests of the nonlocality by means 
of the Bell-type inequality Q] have been made with dif- 
ferent kind of particles including photons Q and mas- 
sive fermions 0- [!■ Q The entanglement is expected 
to play an essential role at quantum phase transitions 
0, where quantum fluctuations manifest themselves at 
all length scales. Several groups investigated this prob- 
lem by stud ying the quantum spin systems (see, e.g., 
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EH 12% Eillh Additional studies have been carried out for 
more complicated systems including both itinerant elec- 
trons and localized spins; the local entanglement for these 
systems have been discussed in a context of the quantum 
phase transitions 0, H3 an d of the Kondo problem (2(| . 
In particular, Anfossi et al |25| performed, within the 
density-matrix renormalization group method, a numer- 
ical comparison between the standard finite-size scaling 
and the local entanglement for the Hubbard model in a 
presence of bond charge interaction (the Hirsch model) at 
the Mott metal-insulator transition. Katsnelson et al |2q ] 
have considered the suppression of the Kondo resonance 
by a probing of the charge state of the magnetic impurity 
which leads to the partial destruction of the entanglement 
between the localized spin and itinerant-electron Fermi 
sea. These examples illustrate a relevance of the concept 
of entanglement for the many-body physics. Moreover, 
the entanglement overwhelmingly comes into play in the 
quantum computation and communication theory |27| . 
being the main physical resource needed for their spe- 
cific tasks. The essential idea is to encode one particular 
qubit, and let it be transported to across the chain to 
recover the code from another qubit some distance away 

HI 

The suppression of the entanglement by decohering ac- 



tions such as noise, measurements, etc., is one of the 
central problems in quantum computation and quantum 
information theory; therefore the concept of entangle- 
ment for mixed states is of primary relevance p9l |. For 
example, it is important to know what happens with the 
quantum computer after the measurement of one qubit 
state; for the case of the quantum system with broken 
continuous symmetry such as Bose-Einstein condensate 
(BEC) or easy-plane antiferromagnet the local measure- 
ments lead to the formation of the "decoherence wave" 
[3fl l3lj . It is interesting to investigate the effect of the 
decoherence wave on the entanglement in the system. 

Motivated by these results, in this paper, we aim on 
the evaluation of the pairwise entanglement in the ID 
Ising-XY model with transverse magnetic field at finite 
temperatures. We identify a region where the thermal 
entanglement is non zero while it is zero at zero temper- 
ature, which results from the entanglement of the excited 
states as it will be explained below (section[nJ . We study 
the dynamical response of the system in the non vanish- 
ing entanglement region, which is the useful region for 
quantum information processing, after a projective mea- 
surement on one local spin. We find that, similar to the 
case of the BEC considered earlier [3(| the spin decoher- 
ence wave appears propagating with the velocity propor- 
tional to the interaction strength. We investigate also 
the zero temperature case studying the time-dependent 
correlation functions and we discuss the relation between 
the dynamics of the magnetization and the entanglement 
(section ITTTj) . Our conclusions are given in section llVl 



II. THERMAL ENTANGLEMENT IN THE 
ISING-XY MODEL WITH TRANSVERSE FIELD 

In this section we present the solution of the N-sites 
Ising-XY model with transverse field following the stan- 
dard method [32l l33l| . We proceed with the Hamiltonian 

H = f>[(l + 7KX +1 + (1 - 7KVf +1 ] + at) , (1) 

i=l 
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Fig. 1: (color online) Pairwise entanglement for the nearest 
neighbors in the isotropic XY model with transverse field at 
the equilibrium as function of /3 and A. Cequ is defined by 
Eq.(12) for the thermodynamic equilibrium state. 




Fig. 2: (color online)Same as Fig. [T]for 0.8 < A < 0.9. 



where af are the Pauli operators, obeying the usual com- 
mutation relations [of ,0^] = 2ie abc 5 v of . The Zeeman 
energy in the external magnetic field, as well as the 
Planck constant H, have been set to 1. We assume cyclic 
boundary conditions, i.e., the index i in the sum £Q| runs 
over 1 . . . TV with Sn+i = Si- This Hamiltonian can be 
diagonalized by means of the Jordan- Wigner transforma- 
tion [32tl33| that maps spins to one-dimensional spinless 
fcrmions with creation and annihilation operators Cj and 
cj. The Hamiltonian Eq.(l) in the fermionic operator 
representation is represented by the quadratic form 

H = Q2 4AiJCj + l^Bijc] + H.c.]) +N, (2) 

where A iyi = -1 and Aij+i = -—A = Ai+ij, Bis+i = 
— 5A7 = —Bi + ij, and all other A+j and B^j are zero. 
The quadratic Hamiltonian J2J can be diagonalized by a 



Magnetization 




Fig. 3: Site magnetization as function of the temperature and 
time with m — I = 1 and A = 0.8. 

linear Bogoliubov transformation of the fermionic oper- 
ators, 

Vk = 22(gkiCi + hki4) , (3) 



vl = ^2(9^1 + h klCi ) , (4) 

i 

where the g k i and hki can be chosen to be real. After 
that it takes the diagonal form 

H = J2 A kvhk - \^k , (5) 

k 

where 



Afc = Vil* sin k) 2 + (1 + A cos k) 2 . (6) 

After the diagonalization of the Hamiltonian now we 
can proceed with the evaluation of the thermal pairwise 
entanglement. The pairwise entanglement, as its name 
indicates, measures how two spins separated by a dis- 
tance r are entangled. This measure is to be accom- 
plished by evaluating the pairwise entanglement of the 
two-site density matrix after tracing out all other spins 
in the chain. We evaluate the entanglement of this states 
by using the concurrence which is defined as [2jj 

C = max{Ai - A 2 - A 3 - A 4 , 0} , (7) 

where A's are the square roots of the eigenvalues in de- 
creasing order of the matrix pab iPy ® a yP*AB a v ® a y)> 
where p AB is the corresponding complex conjugation in 
the computational basis {|++), |+— ), | — +), | )}■ As 
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usual, at the thermal equilibrium the system is described 
by the canonical ensemble density matrix 



-m 



(8) 



where Z = Tre~@ H is the partition function of the sys- 
tem. Thus the reduced two-site density matrix, after tak- 
ing into account the symmetries consideration, assumes 
the following form 



-pH 



pij = 



3 



z 



(<r*)(4 



52(44)4 ® 4) ■ 



I + I®a 3 z ) 
(9) 



k=l 



The correlation functions that show up in the density 
matrix Eq.Q are well known |33l | and for the nearest 
neighbor case considering here these correlation functions 
are giving by 

(ai4) = g. 1 , (44) = g 1 , (aX) = ^ 2 -^-i, 



quantum fluctuations likely to dominate the behavior of 
the system and it is the region where the quantum infor- 
mation processing should be studied since it is well know 
that the entanglement is the main resource for quantum 
information. Thus the study of the dynamics of entan- 
glement for this region is relevant. In the next section we 
will study the response of the system after a projective 
measurement in this region. 



III. SPIN DECOHERENCE AFTER A 
PROJECTIVE MEASUREMENT 



Magnetization 
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(cO = -Go ■ 



(10) 



with 



i r 

Gmi = — I dkcos[k(m — + Acosfc) 
t Jo 



tanh(A fc /3/2) 



A, 



— / dk sin[fc(m — i)] sin k 
n Jo 



tanh(A fc /3/2) 
AT 



(11) 
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Thus the concurrence for the case of isotropic XY model, 
7 = 0, reads 



C = max{0,||£i 



1 



(l + eo 2 -^ 2 ) 2 -^ 2 !}, (12) 



and at T = we have Go = 1 and 



Gi = 



X < 1 
A > 1 



Fig. 4: Site magnetization as function of the temperature and 
time with m — I = 1 and A = 0.99. 



As mentioned above we will study the consequences 
of the local projective measurement and analyze the sys- 
tem behavior in terms of the spin decoherence wave. We 
introduce the operators 



The two-site entanglement between the nearest- 
neighbors for the isotropic XY model is shown in Fig^ 
One can see from this figure that there is a region where 
the entanglement increases with the temperature increase 
whereas for A < 1 the entanglement is zero at zero tem- 
perature and it remains zero until a critical temperature 
where the system starts to be entangled. Fig|5| displays 
with more details a relevant region with 0.8 < A < 0.99. 
It is clearly seen in this figure that there is a strong 
enough entanglement in the region 2 < (3 < 20, however, 
outside this region no entanglement can be observed. 
This entanglement transition should be understood as an 
effect of the entanglement of the higher excited states. 
Note that the ground state in this case is unentangled 
with all spins pointed in the same direction. Similar ob- 
servation has been made in Refs0 and In the region 
where the entanglement does not vanish we expect the 



At 



c\ + Ci 



B, 



(13) 



After a selective projective measurement with the projec- 
(P = pt) j which means that the 



1-A,B, 



imP r\-.— 

positive z direction of the local spin I is the measurement 
result (a general measurement will be considered below), 
the mean value at time t for an operator A reads 



(A(t)) 



TrpPA(t)P 
TrPpP 



(14) 



where A(t) — e A(0)e~ . Since we are interested 
in the evaluation of the time-dependent average value of 
the magnetization in the z direction at site m we have 
A(0) = <j z /2 = -AmBrn, i n order to evaluate (A(t)) 
we write the operator Ai, and B>i in term of r\ and r]\ 
operators using the inverse transformation of Eqs.lO,l@J|. 
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Fig. 5: Site magnetization as function of the temperature 
and time with m — I = 1 and A = 0.8 for the case without 
knowledge of the measurement outcome. 




Fig. 6: Site magnetization as function of the site location 
x = m — I and time at fixed /3 = 10 and A = 2. 

Since the Hamiltonian is diagonal being written in terms 
of r/ operators we have 

4(t)=exp(iA fc t)^(0). (15) 

Thus after a straightforward little algebra we found the 
following expression for A(t) 

A{t ) = T" E a ?i' A i A i' + Pu> A i B i' + lu-BiBv , (16) 

ii' 



Fig. 7: Site magnetization at distance x from the measure- 
ment point. 



where 

®ii' = §miGmi'\ Pa' = ^mi^rni' CmiCmi' j Tzi' = Gmi^rni' i 

(17) 

and in the thermodynamic limit (N — » oo) we have 

1 /" r 

4>y.v = i>nv = - dk cos k(fi - v) cos(A fc t) , (18) 
n Jo 

Gn U = — [ dkcos[k(fi — + Acosfc) — ^ fc ^ 
7t Jo ' " A fe 

_ ^2 r rffcsin[MM _, )]sinfc !i^. (19) 

We define the operator A = J2 U , (3™,AiBi' which corre- 
sponds to the part of A(t) with real coefficients. Clearly 
we have (A)p = (A}p, therefore after some calculations 
we obtain 

= 4 ( g „ 1 + l) + S P«' UiBi^B^p 

ii' 

+ {A i B i ,A l B l ) -{A l B l A i B i ,A l B l ) p }), (20) 

where 

(AiBiAiBi')p = Gh'Gu + SiiSt'i - GuGw , (21) 



{AiB v AiBi) p = Gu'Gu + 5ii5i>i - GuGi'i , (22) 
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{A l B l A l B l , AiBi) p = -45 H ,S u Gu + 8 u {G n + Q w ) + 6 w (G a + Gu) + QuQn{Qu - Gu) + GuGiAGu - Go) - Gu> ■ (23) 

I 



Here we have used the Fermi distribution function the isotropic XY Model, 7 = 0; the case 7^0 will be 

addressed in the future. Thus, the magnetization at the 



e -pA k + i ■ 

For simplicity, in this paper we will further consider only 



(24) site to in the z direction can be written as follow 



J 



(Ah 



-r {G mm + (2Gu + 1)[2(^ - G 2 ml ) + a-a']+ 4(G ml (3' ml - c^ ml f3 ml ) + 2{(3'* l - (3 2 ml )} 



(25) 



where the expressions of a, (3 m h a 'j P' m i are given in the 
Appendix. 

The magnetizations at the neighboring site m — I = 1 
for the cases A = 0.8 and A = 0.99 are shown in Figs[3] 
and^J respectively, as functions of the inverse tempera- 
ture (3 and of the time. The magnetization oscillate in 
time with the frequency proportional to A which is a re- 
sult of the propagation of the decoherence wave after the 
local measurement 30] . One can see that the amplitude 
of these oscillations increases in the close vicinity of the 
quantum critical point A = 1. Figs |3| and ^demonstrate 
also that the amplitude increases with the temperature 
increase. It is connected probably with the thermal en- 
tanglement in the system under consideration (see the 
previous section). 

For a complete von Neumann measurement |27| with- 
out knowledge of the measurement outcome the mean 
value for the operator A at time t is 

(A(t)) = Tr P PA{t)P + Trp(l - P)A(t)(l - P) , (26) 



Thus we have 



(A) 



1 



{Gn 



AGu 



il - G ml) + (« - ot') 



^{G m i(3' ml — 4> m lPml)} ; 



(27) 



The magnetization at the neighboring site m — I = 1 
after the measurement is plotted in Fig [5] for A = 0.8 as a 
function of f3 and the time. One can clearly see from this 
figure that there is a reduction of the amplitude of the os- 
cillation with respect to the selective measurement; this 
effect is due to the mixing of the state as here we don't 
know the results of the measurement outcome. Fig|S]dis- 
plays the propagation of the decoherence wave, that is, 
the magnetization distribution as a function of the dis- 
tance from the measured site x and time t at a fixed value 
of = 10 and A = 2. 

Now we consider the effect of the local measurement 
on the pairwise entanglement in the system. This require 




Fig. 8: Single-site entanglement S at distance x from the 
measurement point. 



the evaluation of a correlations functions such as {afa^). 
Here we will present the results for an interesting partic- 
ular case, namely, for A < 1 and at zero temperature. 
Then, there is no entanglement in the system before the 
measurement since all spin are pointing in the same z- 
direction in the ground state. Therefore a projective mea- 
surement of the z-component of the magnetization will 
not provide us any nontrivial information and will not 
generate entanglement. However, we will show that the 
projective measurement of the x-component does create 
the entanglement. After the projective measurement in 
the x direction at site m with positive outcome the wave 
function will be 



l*)r 



1 



4 



V2 



\vac) . 



(28) 



Fig. 9: Two-site entanglement C between the site m and the 
site i at zero temperature with x = i — m. 



Fig. 11: The fidelity of the quantum channel is shown at site 
i. 



be written in the form 




Fig. 10: (color online) (a z ) (dotted line) and single site en- 
tanglement (solid line) are shown at the site m. 



At time t we have 

4 n {t) = Y J i. G ml + 4>ml 
I 



><4 = X>i(*W 



(29) 



where in the thermodynamics limit wi(t) — J m _/(At) and 
J n (x) is the Bessel function of order n. Thus the time- 
dependent wave function after the measurement will be 



IT tiW 1 +l2l W l(' t ) C l I \ 

l*m(*)) = ^2 |WQC) 



(30) 



Thus, the time-dependent two-spin density matrix can 



Pij 







WiWj 



WiWj 
2 





Wj 

2 



It) 



and the corresponding one-particle density matrix is 
/ w 2 

Pi = 5 ( o ""*„2 ) • The magnetization of the site 



Wi 

i is therefore given by (c a )/2 = Trpcr^/2 = (wf — l)/2. 
The single-site entropy, 5(/0j) = — Trp^ logp^, which char- 
acterizes the entanglement of one spin with the rest of the 
chain, can also be evaluated in this case. The pairwise 
entanglement for the two-site density matrix can be eval- 
uated using Eq. (7J| . A straightforward algebra lead to the 
following expression for the concurrence: 



C 



WiWj . 



(31) 



In Figs[7||BllIIlwe show the single-spin quantum entropy 
and the magnetization at distance x from the measure- 
ment point, as well as the pairwise entanglement between 
site to and i as functions of time and of x = i — to. We 
conclude from these figures that the single-site entan- 
glement and the pairwise entanglement propagate with 
the velocity proportional to the interaction strength, like 
the spin decoherence wave. In Fig llUl we display the site 
magnetization together with the single-site entanglement 
which demonstrates clearly that these two quantities os- 
cillate coherently. This confirms that the dynamics of 
the spin decoherence reflects in some sense the dynamics 
of the entanglement in the system. 

The fidelity of the communication at site m through 
the channel is the probability that a channel output pass 
a test for being the same as the input conducted by some- 
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one who knows what the input was. It can be defined as 



(4>m\Pi\4>r, 



(32) 



where |i/>„ 



|Q)+|i) 

V2 



is the state of the site m right 



after the measurement. In fact, the spin chain acts as an 
amplitude damping quantum channel where the initial 
state is transformed under the action of the superopera- 
tor $ to E3 



p -> %{p) = MopMl + M lP M\ 



with the Kraus operators such Mq 



(33) 



and 



Mi 







o 

7 o 



where as usual Mi describe the 



quantum jump and Mq represent no quantum jump. The 
fidelity of the channel is shown in Fig^J One can see 
clearly from this figure that the channel can be efficiently 
used to transmit the quantum information. The fidelity 
has a maximum value for x — i — m ~ Xt. This means 
that the quantum state is transported with the velocity 
proportional to the interaction strength A similar to the 
decoherence wave. After a time t = x/X the state can 
be recovered with maximum fidelity at a distance x from 
the initial site m. 



one-dimensional Ising-XY model with transverse mag- 
netic field. Our findings indicate that the behavior of 
entanglement with respect to temperature, at least for 
moderate values of temperature, is quite complex. In 
particular, we have found that for some ranges of temper- 
ature, entanglement in the system can grow with increas- 
ing temperature, which results from the entanglement of 
the excited states. We have studied the dynamical re- 
sponse of the system in a relevant region for quantum 
information processing after a projective measurement 
on one local spin which leads to the appearance of the 
"decoherence wave" that propagate with velocity propor- 
tional to the coupling constant A, similar to the case of 
the Bose-Einstein condensate studied earlier 30] . One 
motivation behind our study is to know what happens 
with the quantum computer after the measurement. We 
have investigated for specific case (T = and A < 1) the 
dynamics of the entanglement and the spin decoherence 
wave and we have found that those quantities propagates 
coherently through the chain with the same velocity that 
is proportional to A. The fidelity of the channel has been 
shown to be represented as amplitude damping channel. 
Finally, a generalization to 7 7^ and study of the en- 
tanglement dynamics in such system are desirable. 
Appendix 



IV. CONCLUSIONS 



In this paper, we have evaluated the equilibrium pair- 
wise entanglement at finite temperatures in the isotropic 



Using the standard properties of the Fourier transfor- 
mation, one has 



J 



,ik{i—i ) 



(34) 



1 1 
Putting 7 = wc find 



„ 1 C „ / , N tanh(A fc /3/2) 2/ , . 
a = y j <pmi4>mi>yii> — — / dk(l + Acosfc) cos (Ajfet) 



E 



(35) 



1 r 

bmiGii — — I dkcos[k(m — Z)](l + Acosfc) 

I" JO 



tanh(A fe /3/2) 
A fc 



cos(Afei) 



(36) 



X = / GmiGmi'Gii' = / dk(l 

^7 n Jo 



+ Acosfc) 



3 tanh(A fe /3/2) 2 

- sin (AfciJ 



Af 



(37) 
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P' m i = E G ^i, = - r dk cos[fc(m - 0] (1 + A cos fc) 2 tanh ( A ^/ 2 ) sin(Afct) 
n Jo A k 



k 

(38) 
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